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Abstract. The role of ion-ion correlations in suppressing neutrino-nucleus elastic scattering during stellar core collapse is 
reinvestigated, using two different equations of state. We test the improved description by Itoh et al. against the treatment 
suggested by Horowitz and find that the stronger cross section reduction for small momentum transfer in the former case does 
not lead to noticeable changes of the core deleptonization and entropy increase during collapse, because the improvements are 
relevant below neutrino trapping conditions only for very low neutrino energies, corresponding to a very small phase space 
volume. Treating screening effects for ionic mixtures by the linear mixing rule, applied to the collection of representative heavy 
nucleus, a particles, and free nucleons, which is assumed to characterize the composition in nuclear statistical equilibrium, we 
cannot determine mentionable differences during stellar collapse, because a particles are not sufficiently abundant and their 
coherent scattering opacity is too small. 
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1. Introduction 

Heavy (iron-group or more massive) nuclei dominate the com- 
position in stellar iron cores until nuclear densities are reached 
in the inner core and the bounce shock raises the entropies in 
the outer core to values where free nucleons are favored in 
nuclear statistical equilibrium (NSE). During the infall phase 
therefore coherent, isoenergetic scattering off nuclei is the 
main source of opacity for the elec tron neutrinos p roduced by 
electron captures (cf., e.g.. lBruennlll985l Il989albl) . Neutrino- 
nucleus scattering thus hampers the free escape of neutri- 
nos, is responsible for neutrino trapping around a density of 
10 12 gcrrT 3 , and regulates the deleptonization and increase of 
entropy during core collapse. 

In the medium of the supernova core nuclei are coupled 
strongly with each other by Coulomb forces. They thus form 
a highly correlated plasma, in which the interactions of neutri- 
nos with wavelengths larger than the average ion-ion separa- 
tion aj on (corresponding to neutrino energies e v & 2hc/a\ n „ ~ 
20 MeV) are reduced by phase interference effects llltohl l975). 
Th e correspo nding ion s creening was mor e recently calculated 
bv lHorowitd Jl997l) and lltoh et alJ J2004I). and investigated in 
its effe cts on stellar core collapse bv iBruenn & Me zzacanna 
lll997l) . The latter authors employed the correction factor for 
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neutrino-nucleus scattering cross sections as given by Horowitz 
Jl997l>. 

Iltoh et alJ fcOftfl) . however, pointed out that the Monte 
Carlo calculations, which Horowitz's fit was based on, did not 
allow him to accurately represent the cross section reduction 
for low neutrino energies, i.e., for energies e v % %cla 10B ~ 
10 MeV, thus undere stimating the importance of ion-ion cor- 
relation effects. Ito h et alJl| 2004) provided a more accurate an- 
alytic fitting formula by using the correct behavior of the liq- 
uid structure factor for small momentum transfer in neutrino- 
nucleus scattering. 

The investigations presented in this work have two goals. 
On the one hand we aim at studying the differences for stellar 
core collapse and the formation of the supernova shock which 
arise fro m the improved d escription of ion screening as sug- 
gested by Itoh et al. (2004), compared to a treatment using the 
formulae of iHorowitzl (119971). On the other hand we intend to 
explore the sensitivity of the evolution to ion screening effects 
associated with the ionic mixture of nuclei and nucleons that 
are present during core collapse. In accordance with the treat- 
ment of NSE in current equations of state (EoSs) for super- 
nova simulations, we consider the nuclear components to be 
free neutrons, free protons, a particles, and one kind of heavy 
nucleus which is considered as representative of the NSE distri- 
bution of nuclei beyond 4 He. Two different nuclear EoSs with 
largely different a mass fractions during core collaps e are em- 
ployed . The first EoS ("L&S"), provided by Latt imer & Swestvl 
Jl99ll) . is based on a compressible liquid drop mode l and 
uses a Skyrme force for the nucleon interaction JLattimer et alJ 
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Fig. 1. The ratio of the angle-averaged ion-io n correlation fac- 
tor as calcul ated with the fitting formula o f IHorowitzl dl997h 
(see Fig. 9 in lBruenn & Mezzacappalll997l) relative to t he one 
obtained from the prescription provided bv lltohet all 12004) 
(see their Fig. 2) as a function of the variable £ for different 
values of the dimensionless parameter F. 



Our choice of the compressibility modulus of bulk nu- 
clear matter is 180MeV, and the symmetry energy parameter 
29.3 Me V, but the differences in the supernova evolution caused 
by other values of the c ompressib ility were shown to be minor 
(Thomps on et al J 12003: Swe stv et al.lll994l) . The second EoS 
used here ( "Shen") is the new relativistic mean field EoS of 
Shen et al 1 dl998albl) with a compressibility of nuclear matter 
of 28 1 MeV and a symmetry energy of 36.9 Me V. 

The paper is organised as follows: In Sect. El th e prescrip- 
tions f or t he ion-ion correla tion factor provided by IHorowitzl 
dl997l) and lltoh et al.l (120041) are briefly summarized, in Sect. [3] 
the input (code, initial stellar model) in our simulations is de- 
scribed, in Sect. |4]the results are presented, and in Sect. |5] we 
finish with conclusions. 



2. Ion-ion correlation factor 

In this study we make use of two different fitting formu- 
lae for the angle-averaged correlation factor (Si oa ), which de- 
scribes the reduction of the neutral-current scattering of neu- 
trinos off nuclei by ion-ion correlation effects. It is used as a 
multiplicative correctio n to the neutrino-nucleus iso energetic 
scattering opacity (cf.lHorowitzl ll997: Bru enn & Mezzaca ppa 
1997URampp & .Tankall2002F " The first formula is provided by 
Horowitz! d 19971) and is based on Mon te Carlo results. The sec- 
ond one is given by lltoh et all d2004l) and was obtained from 
data calculate d with the improved hypernetted-chain method 
(see lltoh et al 1 ll983l and references therein) for a classical 
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Table 1. Overview of neutrino-matter and neutrino-neutrino in- 
teractions included in our simulations. For each process we 
provide reference(s) where more information can be found 
about physics and approximations employed in the rate calcu- 
lations. The numerical implementation is descr ibed in detail in 
iRampp & Jankal d2002h and lBuras etail d2003h . The symbol v 
represents any of the neutrinos v e , v e , v^, v^, v T , v T , the symbols 
e~, e + , n, p and A denote electrons, positrons, free neutrons and 
protons, and heavy nuclei, respectively. The symbol N means 
neutrons or protons. 



we also take into acco unt (the rather small) electron screen- 
ing effects according to IHorowitzl \ 1997b by applying the addi- 
tional correction factor of his Eq. (19) to the rates of coherent 
neutrino-nucleus scattering. 



2. 1 . One-component plasma 

If the stellar plasma consists of only one nuclear species of ions 
(Z; on , A; on ), the ion sphere radius which gives the mean interion 
distance is defined as 



Ami-,, 



1/3 



(1) 



where n mn = ri\,Xi on /Ai on is the number density of the ions with 
mass number A; on , charge Z; on e, and mass fraction X lon («b is 
the number density of baryons). The strength of the ion-ion 
correlations is characterized by the dimensionless parameter 
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one-component plasma. In all simulations with ion screening 



which is the ratio of the unshielded electrostatic potential en- 
ergy between the neighbouring ions to the thermal energy. 
Here T\o denotes the temperature in units of 10 10 K and p\i 
the mass dens i ty in 10 12 g cm' 3 . Note that the definition of 
r in IHor owitzl fl 997l) differs from the ones used h ere and in 
Bruenn & Mezzacappa ill 9971) and lltoh etal J §004) by a fac- 
tor (4tt) in the denominator of Eq. (TJ}- In Horowitz's notation 
this factor is absorbed in the employed value of e 2 . 
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Fig. 2. Composition as a function of density at the center (top) and at an enclosed mass of 0.5 M (bottom) for the evolution during 
core collapse with the EoSs of Lattimer & Swesty (1991; left) and Shen et al. (1998a,b; right). The plots show the evolution of 
the number fractions Y for free neutrons, protons, a-particles, and the representative heavy nucleus until shock formation. The 
results were obtained with ion-ion screening treated according to case "EoSjonJtoh" (see text). 



iHorowitj i 19971) provided the following fitting formula for 
the angle-averaged suppression factor (Si on (g, r)) Hor : 



(S ion(£, 0) H or 



and 



( 6 



1 + exp 



V i=0 



for 3 + ^,(3) 



<Sion(£,r)>Hor = 1 - otherwise , 



(4) 



to be app lied for 1 < T < 150; for values of T < 1 or 
F > 150. iHorowit d l l 19971) recommends to simply set T to 1 
or 150, respectively. The /?, are coefficient functions of Y deter- 
mined from fits to Monte Carlo data. In Eq. the variable £ is 
the ratio of the mean ion-ion separation a lon to the wavelength 
for neutrinos (during core collapse primarily electron neutrinos 
with energy e v ), i.e., 



He 



(5) 



lltoh et alJ d2004 provide a different fitting formula, see 
Eqs. (23)-(26) in their paper, which is restricted to the case of a 
strongly degenerate electron gas, a usually well fulfilled condi- 
tion during stellar core collapse. Their treatment gives different 
results for the ion-ion correlation factor, (S i n(£> OX m the limit 
of low neutrino energies (£ $ 1). This can be seen in Fig.^ 
lltoh et all J2004I) argue that the reason for t his difference com- 
pared to the description by Horowitz ( 1997) is their correct cal- 
culation of the liquid structure factor S (k) for small momentum 
transfer k. This makes the suppression of neutrino-nucleus scat- 
terin g by ion-ion corr elations more important than es timated by 
lHorowitzNl997l) and lBruenn & Mezzacappal i l 19971) . 

The simulations presented here intend to study the dynam- 
ical consequences of these differences during supernova core 
collapse. 

Up to now we have considered a stellar medium consisting 
of only one nuclear species and referred to formulae derived for 
a classical one-component plasma. However, in the collapsing 
core of a massive star a mixture of nuclei besides free neu- 
trons and protons is present, usually approximated by a parti- 
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cles plus one representative, heavy nucleus. To deal with that 
we decided, in the simplest approach, to calculate the suppres- 
sion factor (S j on (£, T)) for a particles and for the heavy nucleus 
independently, using the values for the average distance a\ on be- 
tween ions of the same kind, derived from the number densities 
n a or ha, respectively. This assumes that different ionic com- 
ponents coexist without collectively affecting the screening of 
neutrino-nucleus interactions. Alpha particles thus change the 
ion screening for heavy nuclei only by the fact that their pres- 
ence may reduce the number density of heavier nuclei. 

2.2. Ionic mixtures 

For a liq uid mixture of different ions (Zj,Aj) including free 
protons, fitoh et alJ ( 120041) suggest a modi fied treatment, refer- 
ring to earlier work by Itoh et al Employing the so- 
called linear mixing rule, one can extend the calculations of 
neutrino-nucleus scattering cross sections, obtained for a one- 
component ion liquid, to the case of multi-component fluids. 
The ion sphere radius for an ion j in the mix is now given 

by 



1 J 



1/3 



with a e being the electron sphere radius, 



1/3 



1,471- Z/«< 



(6) 



(7) 



where the sum extends over free protons and all nuclei with 
number densities n,. The dimensionless variable then be- 
comes 



(8) 



The ion-ion correlations of nuclear species j depend on the di- 
mensionless parameter Fy defined as 
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The angle-averaged ion-ion correlation factor (S i on (^j,Tj)) is 
now evaluated for a particles and hea vy nuclei with the fitting 
formula provided bv lltoh et alJ J2004). using fjj and Tj as given 
in Eqs. (|8j and (j9jl, respectively. 

Applying the naive procedure of Sect. l2.1l for the conditions 
in a supernova core shows that usually £ » 1 for a particles be- 
cause a's are less abundant than heavy nuclei in the central part 
of the core and during most phases of the collapse. Their T is 
then less than unity. Both factors diminish ion screening for a's 
to a negligible level. In contrast, following the description in 
this section, the presence of a's can affect also the ion screen- 
ing of heavy nuclei by reducing the interion separation (Eq.|6j 
and thus (Eq. |8j. Moreover, Tj for heavy nuclei according 
to Eq. (|9} might become larger than in Eq. (0. Therefore the 
presence of a's has an indirect influence on neutrino scattering 
off heavy nuclei and thus on ion-ion correlations during stellar 
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Fig. 3. Central electron fraction Y e and lepton fraction Y\ ep as 
functions of central de nsity for core collapse calculations with 
the EoS of lShen et al1dl998albl) . Model Shen _ion_off does not 
include the effects of ion-ion correlations, Model Shen jon_Hor 
uses t he description of ion screening according to iHorowitzl 
lll997l) . Model Shen ion Jtoh employs Itoh et al.'s (2004) treat- 
ment for a one-component plasma, and Model ShenJon_mix 
their treatment of ionic mixtures. 



core collapse, despite the fact that the screening effects for a 
particles a re still sm all because their T is usually below unity 
(following Horowitz ( 1997), F is then set to unity for evaluating 
the angle-averaged cross section suppression factor). 

While our procedure for treating the effects of a parti- 
cles in an ionic mixtur e with heavy nuclei and free nu cle- 
ons adopts the recipe of lltoh et alJ ( 12004 . ISawveJ J2005I) re- 
cently discussed an alternative approach to the problem for 
multi-component fluids by applying the Debye-Hiickel approx- 
imation. He pointed out that in multi-component plasmas the 
ion-ion correlation effects might be greatly reduced, leading 
to much larger neutrino opacities than for a one-component 
plasma, even if the constituent ions have only a small range of 
N/Z ratios. Electron density fluctuations for an ionic mixture 
enhance this tendency. Sawyer's calculations therefore yield a 
result which is opposite to our application of the linear mixing 
rule for calculating multi-component plasma parameters. They 
tend to bring one back closer to the case with ion-ion corre- 
lations being ignored, simulations of which will be presented 
below, too. Thus we provide a set of models with the intention 
to encompass and bracket the "extreme" possibilities discussed 
in the literature. 



3. Numerical simulations 

The core collapse simulations performed in this work were 
carried out with the neutrino-hydrodynamics code Vertex 
in spherical symmetry . The code is described in detail in 
iRampp & Jankal J2002t) . 

Vertex employs the Prometheus hydrodynamics code 
jFrvxelletalJl 989). which is a direct Eulerian implementation 
of the Piecewise Parabolic Method. It uses an explicit Godunov 
type scheme with second order temporal and third order spatial 
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Fig. 4. Same as Fig.[3]but for the central (matter) entropy s. 
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due to the weak magnetism of nucleons dHorowitzll2002l) . the 
possible quenching of the axial-vector coupling in nuclear mat- 
ter JCarter & Prakashll2002l) . and the reduction of th e effective 
nucleon mass at high densities jReddv et alj [lQ99). Electron 
captures on nu clei are implemented a ccording to the improved 
treatment of Lan ganke et al.l J2003h in regions where NSE 
holds, taking into account the collective e-captures of a large 
sample of nuclei in NSE with rates determined from shell 
model Monte Carlo calculations; the prescription of iBruennl 
(1985) is used in regions which are out of NSE. With this input, 
the production of v/ s by nuclei dominates the one by protons 
during core collapse (Langanke et all2003l) . 

The core collapse simulations presented in this paper were 
started from the 15M Q progenitor sl5a28 from iHeger et alJ 
(2001). 
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Fig. 5. Gas entropy vs. enclosed mass at the moment of shock 
formation in the models of Figs. [3]and0] The shock forma- 
tion is defined by the instant when the entropy in the core first 
reaches a value of 3 k% per nucleon. 

accuracy, which allows for the simulation of high Mach num- 
ber, self-gravitating flows with discontinuities and can handle 
complicated equations of state. The Vertex transport solver 
employs a variable Eddington factor algorithm in spherical co- 
ordinates. Angular moments of the 0(v/c) comoving frame 
transport equation form a system of lower-dimensional mo- 
ments equations, which is closed with a variable Eddington 
factor. The latter is determined from a (simplified) "model" 
Boltzmann equation. General relativistic gravity and corre- 
sponding effects in the transport are taken into account ap- 
proxim atel y. Details about this can be found in Rampp & Janka 
(l2002h and lMareketail i2005h . A comparison with fully rela- 
tivistic calculations revealed excelle nt agreement during core 
collapse and shock formation ( Liebendo erfer et al. 2005). 

The neutrino interactions used in the simulations of the 
present work are summarized in Table [2 Note that neutral- 
current scatterings of neutrinos off nucleons and charged- 
current /3-processes include the effects of nucleon recoil, ther- 
mal motions, an d phase space blocking, nucleon correlations 
in dense media jBurrows & Sawverl ll998, 1999), corrections 



For describing the thermodynamics and composition of 
the stellar plasma, the EoS of | Lattimer & Swestvl dl99ll) and 
the one of IShen et al.l Jl998albl) are applied at high densities 
(p > 6.7 x 10 7 gcirT 3 or p > 2.7 x 10 8 gcirT 3 , respectively). 
At lower densities the EoS contains a mixture of electrons, 
positrons, photons, nucleons and nuclei, with the nuclear com- 
position being described by a simple approximation to a four- 
species NSE for temperatures above about 0.5 Me V. Below that 
temperature the composition is adopted from the progenitor 
star and modified if nuclear burning plays a ro le during col - 
lapse (for details, see Appendix B in iRampp & Jankal l2002). 
The two EoSs show major differences in the abundances of a 
partic les, which can be larg er by up to a factor of ~ 10 in case 
of the Sh en et alJ (11998a b) EoS. This is visible in Fig.[2]where 
the number fractions of free neutrons, protons, ar's and of the 
representative heavy nucleus (whose mass and charge numbers 
typically grow with density until nuclei disappear at the phase 
transition to nuclear matter) are displayed as functions of in- 
creasing density during collapse both at the stellar center and 
at an enclosed mass of 0.5 M Q . Although th eir mass fraction is 
much lower, a particles in the lShenetal"llll998albl) EoS can 
become equally or even more abundant (by a factor up to about 
two) than heavy nuclei in the outer layers of the collapsing core, 
in particular exterior to M(r) ^ 0.5-0.6 Mq. 

For each of the employed EoSs four core collapse simu- 
lations were performed, all starting from the onset of gravi- 
tational instability and carried on until the moment of shock 
formation. Calculations with ion-ion correlation (and elec- 
tron screening) effects in neutrino-nucleus scattering being 
switched off, i.e. for (5; on ) = 1 (and/? e s 1 instead of Eq. 19 of 
Horowitz 1997), are denoted with "EoS_ion_off", where EoS 
stands for "L&S" or "Shen". They are compared with simu- 
lations (models "EoSion_Hor") where the ion-ion correlation 
factor of Horowitz (1997; Eqs. l3l4> is used, and with models in 
which ion-ion correlations are described according to lltoh etall 
J2004I) (models "EoS Jon Jtoh"). Finally, the sensitivity of stel- 
lar core collapse to the treatment of ion screening for ionic mix- 
tures is investigated by simulations (models "EoSJon_mix") in 
which the correction factors (S i on (^j,T ;)) are calculated from 
Itoh et al.'s (2004) formulae with and T ; as given in Sect. 12. 21 
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Fig. 6. Spectra of neutrino energy density (solid line) and energy flux (dashed) for Model Shen_ion_off (left) and 
Model Shen Jon_Hor (right) when a density of 10 12 g crrT 3 is reached at an enclosed mass of 0.3 M Q . 
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Fig. 7. Profiles of electron fraction Y e , lepton fraction 7i ep (left), and (gas) entropy s (right) versus enclosed mass for 
Models Shenjon_off (dotted), Shen Jon _Hor (solid), ShenjonJtoh (dashed), and Shen_ion_mix (dash-dotted) at the time when 
the central density has reached a value of 10 14 g cm -3 . 



4. Results 

Figures [3] and 0] show electron fraction Y e , lepton frac- 
tion Kjep and (gas) entropy s, respectivel y, at the core cen- 
ter during collapse simulations with the IShen et al.l ( 1998a b) 
EoS. For both EoSs employed in this work, the same rela- 
tive changes are found when models without ion screening 
are compa red with ca lculations with ion-ion correlations ac- 
cording to iHorowita (119971) (see Models Shen JonJHor and 

Shen Jon_off in Figs.RlandRJ. We shall mos tly concentrate here 

— I 1 r i-i 

on the results obtained with the Shen et al. ( 1998a b) EoS, be- 
cause a particles are much more abundant there (see Fig. |2) 
and m any aspects of ion screening in simulations with the 
EoS of lLattimer & Swestvl ill 99 lb were already discussed by 
iBruenn & Mez zacappa ( 1 19971) . Our results agree qualitatively 
with tho se of the latter paper. Quanti tative differences com- 
pared to Bruenn & Mezzacappa (1997) are caused by the in- 
clusion of improved electron capture rates on nuclei in our 



work, which significantly increase electron captures above a 
few 10 10 gcm~ 3 so that lower values of 7ip n and s result af- 
ter tr apping (cf. Lang anke et aljE003t iMartmez-Pinedo et alJ 
120051) . 

As explained in detail by IBruenn & Mezzac appa dl997l) . 
the screened cross section for neutrino-nucleus scattering re- 
duces the transpor t optical depth of low -energy neutrinos (cf. 
Fig. 3 in Bruenn & Mezzacappa! 1 19971) and allows them to 
escape from the core more easily. This is obvious from a 
flux enhancement of neutrinos at energies e v S5 10 MeV in 
Fig. |6] where the situation is displayed at a density of p = 
10 12 gcm~ 3 . Ion-ion correlations thus cause a decrease of Y e 
and Fjep that is stronger by about 0.02 until neutrino trapping 
sets in (Fig.0. The homologously collapsing stellar core corre- 
spondingly shrinks and the shock forms at a somewhat smaller 
enclosed mass (Fig. [5}. The shock formation is defined by the 
moment when the postshock entropy first reaches a value of 
3 &b per nucleon. A part of the effect visible in Fig. |5] might 
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Fig. 8. The energy source term for neutr ino-electron scatterin g for densities of (a) 10 11 gem 3 and (b) 10 12 gcirT 3 at the stellar 
center in the collapse models with the IShen et alJ Jl998albh EoS. Negative values mean that neutrino energy is "absorbed" 
(net scattering out of the corresponding energy bin), positive values mean "emission" of neutrino energy (i.e., net scattering of 
neutrinos into the energy bin). 



therefore be a consequence of the slightly higher core entropy 
after neutrino trapping in models with ion screening (see Figs.0] 
and|7}. 

Neutrino-electron scattering is very efficient in downscat- 
tering neutrinos from the high energies, where they are mostly 
created by electron captures, to lower energy states. Therefore 
the phase space at low energies is quickly refilled. Figure [5] 
shows the source term for energy redistribution by neutrino 
scatterings off electrons for two density values below trap- 
ping conditions. The downscattering of high-energy neutrinos 
explains why the local energy spectra, dE v /de v with E v be- 
ing the neutrino energy density, are essentially the same in 
Models Shen Jon_off and Shen ionJHor, despite of clear dif- 
ferences between the energy flux spectra of both runs (Fig.|6j. 

On their way out escaping neutrinos transfer a part of 
their energy to electrons in collisions, thus heating the stel- 
lar medium (Bruenn 1986). Since ion screening of neutrino- 
nucleus scatterings reduces the transport opacity and therefore 
the effective optical depth for energy exchange with the stel- 
lar background mainly for low-energy neutrinos, but hardly 
changes the downscattering probability of high-energy neu- 
trinos (Fig. |9j, the larger loss of lepton number leads to an 
increase of the central entropy by about 0.12£b per nucleon 
(Fig-EJ- In Fig.|5]the optical depth t at the center is calculated 



tie) 



= £ 

Jo 



dr (A eS (e, r)) 



(10) 



where A e g (e, r) is the effective mean path for energy exchange, 
i.e., the average displacement between two reactions with en- 
ergy exchange between n eutrinos and target particles (see 
Rvbicki & Lightm anll97l . It can be expressed as 



A e ff(e, r) = ^A t (e,r) A e (e,r) , 



(11) 



when A e is the mean free path for reactions with energy ex- 
change, i.e. neutrino absorption and neutrino-electron scatter- 
ing in the present context, and A t is the total mean free path 



for momentum transfer ("transport mean free path"), which in- 
cludes all processes by which neutrinos interact with the stellar 
fluid. 

Figure ^| displays the ion screening factors (Sj n) super- 
imposed to the normalized neutrino energy spectra dE v /de v 
at the stellar center for densities of 10 11 gem" 3 , 10 12 gcm~ 3 , 
10 13 gcirT 3 , and 10 I4 gcirT 3 . Figure ITUl once more demon- 
strates that the reduction of neutrino-nucleus scattering mostly 
affects neutrinos at energies below the spectral maximum for 
all plotted cases. In combination with Fig. |9] it also shows 
that the trapping conditions for the bulk of the neutrino spec- 
trum are not influ enced strongly by ion sc reening. This was 
identified by Bruenn & Mezzacarjpa ( 1997) as the reason why 
ion-ion correlations have no dramatic effect on the core delep- 
tonization. 

From Fig. [2 it is clear that differences betw een the treat- 
ments of ion-ion correlations by Horowitz ( 1997|) and 
(2004) are largest for £ £ 0.5. The improvements by 



Itoh et al 



Itoh et al 



(2004) are therefore most important for the lowest neutrino en- 
ergies in the energy window affected by ion screening. Since 
the phase space available at such low energies is small, one 
cannot expect large quantitative consequences for stellar core 
collapse. This is confirmed by Figs.[3H5]and Fig.Q Itoh et al.'s 
(2004) description (in Model Shen ion Jtoh) leads to values of 
Y e , Tiep, and s after trapping which are essentially indistinguish- 
able from those obtained with Horowitz's (1997) formulae, 
consistent with the insignificant differences between Models 
Shen _ion_Hor and ShenjonJtoh seen in the other plots. Note 
that the crossing of the Y e -, Ti ep - and s-profiles for simulations 
with and without ion s creening at 0.45 M n (F ig. |7} was also 
present in the results of Bruenn & Mezzacappa ( 1997). Finally 
testing the sensitivity of the core collapse evolution to the treat- 
ment of ion screening for the ionic mixture of free protons, a 
particles and a representative heavy nucleus, we also could not 
discover any differences of relevance. 
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Fig. 9. The optical depth for energy exchange between neutrinos and stellar plasma as a function of the neutrino energy at the 
center of the iron core for densities of (a) 10 11 gem" 3 , (b) 10 12 g cm" 3 , (c) 10 13 g cm" 3 , and (d) 10 14 g cm" 3 . The results are taken 
from collapse calculations with the Sh en et al.l 41998a b) EoS. The left panels show enlargements of the low-energy window 
where ion-ion correlations have the largest effect. 



A more detailed analysis reveals the reasons for this insen- 
sitivity, which are valid for both employed EoSs: Below the 
neutrino trapping regime (i.e., for p 5= 10 12 gcm" 3 ) even for 
low-energy neutrinos (e v * 5 MeV) the parameter £ is larger 
than or around unity (except for neutrinos interacting with a 
particles in an ionic mix when % a is computed from Eqs.|6j-[HJ- 
Moreover, F 50 holds at the same time, implying that the 
ion-ion correlation factors (S\ on {^, F))Hor and (S\ on {^, r))i to h for 
neutrino scattering off heavy nuclei are essentially the same 
(see Fig. 0. Only at densities above the trapping density the 
value of £ drops significantly below unity and F exceeds 50, 
causing visible (typically factors 2-3 for e v % 5 MeV; Fig. llOb 
differences in the ion-ion suppression factors (5; on ). At these 
densities, however, the exact value of the neutrino-nucleus scat- 
tering cross section has no noticeable influence on the evolution 
of the core properties and on the neutrino transport. 



Alpha particles in the inner core do not become sufficiently 
abundant to cause mentionable differences. Their indirect ef- 
fect on ion-ion correlations of heavy nuclei in an ion mixture 
by reducing aj (and thus gj) and increasing Fj (Sect. l2.2t for the 
heavier nuclei is essentially negligible, because their contribu- 
tion to the sum in Eq. is diminished by their number density 
being multiplied with a factor Z a /Zj <K 1. Moreover, a parti- 
cles do not account for a significant contribution to the total 
neutral-current scattering opacity, because the opacity for co- 
herent scattering of neutrinos by nuclei (Z, Af, A) scales roughly 
with N 2 /A and therefore is much smaller for a particles than 
for heavy nuclei. For this reason the direct influence of a par- 
ticles and thus of the suppression of their coherent (elastic) 
scattering cross section for low neutrino energies is miniscule, 
despite the fact that £, a drops below unity already at densities 
p % 10 11 g cm" 3 in a mixture with heavy nuclei. (On the other 
hand, T a turns out to be always less than unity and, following 
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Fig. 10. The cross section suppression factor (S l0 n) for neutrino scattering off heavy nuclei (bold lines) as function of neutrino 
energy at the center of the iron core for densities of (a) 10 11 gem" 3 , (bj 1 12 gcirT 3 , (c) 10 13 gcirT 3 , and (d) 10 I4 gcirT 3 . The 
results are taken from collapse simulations with the IShen et alJ Jl99 8a b) EoS. Also plotted (thin solid lines with crosses) are 
the local spectra of the neutrino energy density (normalized to the spectral maximum), which are practically identical for all 
simulations. 



lHorowitzNl997l) . is therefore set to T a — 1 for evaluating the 
angle-averaged ion screening correction factor.) 

5. Conclusions 

In this paper we presented results from simulations of stel- 
lar core collapse with the aim to investigate the consequences 
of ion-ion correlations in neutrino-nucleus scattering, compar- 
ing Ito h et al.'s (2004) improved description with an older 
one by l Horowitzl dl997h . We employed the EoS of IShen et alJ 
Jl998albh in addition to Lattimer and Swesty's (1991) EoS 
and treated electron c aptures on heavy nuclei according to 
lLanganke et al.1 J2003), making nuclei dominant over protons 
in producing neutrinos up to the density of the phase transition 
to nuclear matter. 

Despite these differences in the input physics, our 
models essentially confirmed the previous calculations by 



iBruenn & Mezz ac anp a 1 1997b who followed lHorowitzl J 19971) 
in their description of ion screening. Because ion screening 
is effective only in a low-energy window where the available 
phase space is rather small, the influence of ion-ion correla- 
tions during stellar core coll apse and on the formation of the 
supernova shock is moderate (Bru enn & Mezzacappall997l) . 

We found that the improvement bv lltohet alJ f2004) does 
not lead to any noticeable differences because it affects only 
neutrinos of very low energies (55 MeV) before trapping den- 
sities (p as 10 12 gem -3 ) are reached. Differences at larger neu- 
trino energies occur only at higher densities and thus do not 
affect the deleptonization and entropy evolution. Effects due to 
the ionic mixture of free protons, a particles, and a representa- 
tiv e heavy nucleus — using the linear mixing rule as suggested 
bv lltohet alJl2004 — were found to be negligibly small, too, 
mainly because the abundance of a particles in the inner re- 
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gions of the collapsing stellar core is too low to affect the ion 
screening of heavy nuclei indirectly (see Sect. 12. 21 . Alpha par- 
ticles do not contribute to the total opacity for elastic neutrino- 
nucleus scattering on a level where their ion screening (which 
becomes sizable only when the mixture effects of Sect. l2.2l are 
accounted for) might be relevant. 

Improving the description of ion-ion correlations for the 
complex mix of heavy nuclei with a large variety of compo- 
nents, alpha particles, and free nucleons in the supernova core, 
however, is desirable. Referring to multi-component calcula- 
tions based on the Debye- Hiickel approximation in the limit of 
small momentum transfer, Sawver<2005) argues that a range of 
N/Z ratios in ionic mixtures can protect against the strong ion 
screening suppression of neutrino-nuclei scattering predicted 
by the effective averages of one-component plasma parameters 
applied in the current literature and in this work. Moreover, 
the description of nearly free nucleons and nuclei in NSE is 
expected to hold only up to a density of about 10 13 gcrrT 3 . 
Above this density and below the normal nuclear matter satura- 
tion density, a pasta phase may develop with nucleons clustered 
in subtle and complex shapes. Correlation effects for coherent 
neutrino scattering can then not be treated within the s ingle 
heavy nucleus appro ximation (e.g., iHorowitz et al.l l2004albt 
IWatanabe et al.l2004 . 
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